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Abstract 

In view of new precise measurements of the B D*£ v decay rate near 
zero recoil, we reconsider the theoretical uncertainties in the extraction of 
I Vct\ using heavy quark symmetry. In particular, we combine our previous 

estimate of 1/ mg corrections to the normalization of the hadronic form 
factor at zero recoil with sum rules derived by Shifman et al. to obtain a 
new prediction with less theoretical uncertainty. We also summarize the 
status of the calculation of short-distance corrections, and of the slope of 
the form factor at zero recoil. We find J-{1) = r]A_(,{l) = 0.93 it 0.03 and 

= 0.7 lb 0.2. Combining this with the most recent experimental results, 
we obtain the model-independent value | Vcb\ = 0.040 it 0.003. 
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1 Introduction 



With the discovery of heavy quark symmetry (for a review see Ref. [|T| and ref- 
erences therein), it has become clear that the study of exclusive semileptonic 
B D*i V decays close to zero recoil allows for a reliable determination of the 
CKM matrix element Vch-, which is free, to a large extent, of hadronic uncer- 
tainties 0-0] • Model dependence enters this analysis only at the level of power 
corrections, which are suppressed by a factor of at least (Aqcd/''tT'c)^- These cor- 
rections can be investigated in a systematic way using the heavy quark effective 
theory ||^. They are found to be small, of the order of a few per cent. 

Until recently, this method to determine | Vch\ was limited by large experimen- 
tal uncertainties of about 15-20%, which were much larger than the theoretical 
uncertainties in the analysis of symmetry-breaking corrections. However, three 
collaborations have now presented results of higher precision It is thus 

important to reconsider the status of the theoretical analysis, even more so since 
the original analysis of power corrections in Ref. has become the subject of 
some controversy |10| . 



Besides reviewing some of the existing calculations, the main purpose of this 
note is to propose a "constructive synthesis" of the two approaches that have 
been suggested to obtain an estimate of the power corrections to the decay form 
factor at zero recoil. These corrections are parametrized by a quantity 5i/m2. The 
"exclusive approach" of Falk and myself has the advantage that it provides 
an exact expression for (5i/m2 involving five hadronic parameters, which are de- 
fined in terms of matrix elements of higher- dimensional operators in the heavy 
quark effective theory. The final numerical estimate is model-dependent, since 
four of these five parameters are not precisely known. The "inclusive approach" 
of Shifman tt al. |]10| provides an upper bound for 5i/m2 in terms of only two 
parameters; however, it is not clear to which extent this bound is saturated. 
We shall combine the two approaches and derive non-trivial constraints on the 
hadronic parameters in the formula for 61/^2- These constraints help to reduce 
the theoretical uncertainty. 

Let us start with a short discussion of the decay kinematics [Q] . The hadronic 
matrix element describing the decay process B D*ii? can be parametrized 
by invariant helicity amplitudes corresponding to transverse and longitudinal 
polarization of the D* meson. As kinematic variable, we choose the product of 
the meson velocities, w = vb-vd*, which is related to the momentum transfer 
to the lepton pair by 

w = . (11 

The differential decay rate dT/dw is proportional to the sum over the squared 
helicity amplitudes, which up to a kinematic factor defines the square of a function 
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J-'{w). The resulting expression is 



dr G 
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^ {niB — rriD'Y ■m%*y/w'^ — l{w + 1)^ 

K^r-^'H. (2) 



Aw m\ — 2wmBrnj:i* + m'\),'^ 



w + 1 {niB — mD*y 

The heavy quark effective theory allows the factorization of the short- and long- 
distance contributions to J^{w) into a perturbative coefficient rj^ and a hadronic 
form factor ^{w): 

T{w) = r,Ai{w). (3) 

In the heavy quark limit, this hadronic form factor agrees with the Isgur-Wise 
function ^{w) 0, 0. We use the notation i{w) to indicate that the two functions 
differ by terms suppressed by inverse powers of the heavy quark masses. Luke's 
theorem determines the normalization of S,{w) at zero recoil (w = 1) up to second- 
order power corrections [Q: 



e(l) = 1+51/^2. (4) 

The strategy proposed in Ref. 0] is to measure the product \ Vcb\J^{w) as a 
function of w, and to extrapolate it to w = 1 to extract 

I Kfel ^(1) = I Vcb\ VA (1 + = I Vcb\ {l + 0[asimQ), 1/my } , (5) 

where we use tuq as a generic notation for rric or mh. The task of theorists is to 
provide a reliable calculation of the symmetry-breaking corrections contained in 
7] A and Sl/m^2 in order to turn this measurement into a precise determination of 
\Vcb\- In Sect. H, we briefly review the status of the calculation of short-distance 
corrections. A new theoretical analysis of power corrections is given in Sect. ^ 
In Sect. ^, we give a theoretical prediction for the slope of the form factor ^{w) 
at zero recoil. This parameter is important for the extrapolation of experimental 
data to w = 1. Section |^ contains a summary of the numerical results and some 
conclusions. 



2 Calculation of r}A 

The short-distance coefficient tja takes into account a finite renormalization of 
the axial vector current in the region rrib > fi > rric. Its calculation is a straight- 
forward application of QCD perturbation theory. At the one-loop order, one finds 

as frub + mc , rub 8\ 
r/A = 1 + — In . (6) 

TT Krub — rric rric oj 

The scale of the running coupling constant is not determined at this order. Choos- 
ing as between asijnb) — 0.20 and asirric) — 0.32, and using rric/mb = 0.30±0.05, 
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one obtains values in the range 0.95 < rjA < 0.98. The scale ambiguity leads to 
an uncertainty of order ArjA ~ [(as/vr) ln(m(,/mc)]^ ~ 2%. 

The calculation can be improved by using the renormalization group to re- 
sum the leading and next-to-leading logarithms of the type [a^ ln(m6/mc)]", 
as[as ln(mb/mc)]", and {mc/mb)[as ln(mb/mc)]" to all orders in perturbation the- 
ory ||15|-||18|. A consistent scheme for a next-to-leading-order calculation of rjA 



has been developed in Ref. [0]. The result is 



6/25 /i ^ 1 p^fii a,(m,) -a,(mb) _ 8a, (mj 

vr 377 



rrib V54 27 18 25 

I 2Q;,(m) mg ^^mb| 

where a; = as(mc)/as(mb), and > m > m-c. The numerical result is very 
stable under changes of the input parameters. Using Ky^ = (0.25 ± 0.05) GeV 
(for four flavours) and rric/mb = 0.30 ±0.05, one obtains rjA = 0.985 ±0.006. The 
uncertainty arising from next-to-next-to- leading corrections is of order ArjA ~ 
{as/irf ~ 1%. 

Equation (0) is an exact result to a given order in perturbation theory. We 
stress that, since next-to- leading effects are properly included, it is not only valid 
for large values of ln{mb/mc)- Therefore, we disagree with the criticism of this 
calculation by the authors of Ref. 0]. Of course, it would be desirable to know 
the non-logarithmic terms of order in rjA, but we see no reason why these 
terms should be unusually large. Taking this usual perturbative uncertainty into 
account, we believe it is conservative to increase the error by a factor 2.5 and 
quote 

r/A = 0.985 ± 0.015 . (8) 



3 Anatomy of ^1/^2 

Hadronic uncertainties enter the determination of | Vcb\ at the level of second- 
order power corrections, which are expected to be of order (Aqcd/'^c)^ ~ 3%. 
For a precision measurement of | Vcb\, it is important to understand the structure 
of these corrections in detail. In our discussion (as in all previous analyses), we 
will investigate the 1 / mq corrections at the tree level, thereby neglecting effects of 
order as{mQ)/mQ. In particular, we will not discuss the running of the hadronic 
parameters of the effective theory. In view of the theoretical uncertainty in the 
estimate of these non-perturbative parameters, this is a safe approximation. 
Using the technology of the heavy quark effective theory, Falk and myself have 
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derived the exact expression 

^1/...' = -f^ - ^If:^ - ^) + f^i + - ) . (9) 



\2mc 2mi,J \2mc 2mf,/ Arricmh \3 

which depends upon five hadronic parameters that are independent of the heavy 
quark masses. They have the following physical significance: ip and iy parametrize 
the deficit in the "wave-function overlap" between b- and c-fiavoured pseudoscalar 
(P) and vector ( V) mesons. For instance, ip is defined as 

{D{v)\ c% \B{v)) = 2v%y {l - - ^y^P + 0{l/ml)^ , (10) 

where riv — 1.03 is a short-distance correction factor |l|], and we use a mass- 
independent normalization of meson states. The corresponding relation for vector 
mesons defines iy- The parameter Ai = ~{Pq) is proportional to the kinetic en- 
ergy of the heavy quark inside a heavy meson, and A2 = ^irriy — mp) determines 
the vector-pseudoscalar mass splitting arising from operators in the effective La- 
grangian that break the heavy quark spin symmetry. From the observed mass 
splitting between B and B* mesons, one obtains A2 — 0.12 GeV^. Finally, \q2 
parametrizes certain matrix elements containing two insertions of operators that 
break the spin symmetry. In our analysis below, we will assume that this pa- 
rameter is small, i.e. of a magnitude similar to A2 or smaller. This assumption is 
supported by QCD sum rule calculations of other spin-symmetry-breaking cor- 



rections to heavy quark decay form factors |2^, |21| . 

With the exception of A2, estimates of these hadronic parameters are model- 
dependent. In Ref. 0, we made the simplifying assumptions that ip = iy, and 
that the corrections represented by Xq^ are negligible. Using then reasonable 
values such as ip = iy = (0.35 ± 0.15) GeV^ and -Ai = (0.25 ± 0.20) GeV^, one 
obtains ^1/^2 = —(2.4 ± 1.3)%. Here and in the following, we take rrib = 4.80 
GeV and rric = 1.45 GeV for the heavy quark masses. In Ref. [Q], the error in the 
estimate of Si/m^ has been increased to ±4% in order to account for the model 
dependence and higher-order corrections. A very similar result, —5% < Si/^n^ < 0, 



has been obtained by Mannel ||2^ 



Recently, Shifman et al. have suggested an alternative approach to obtain an 
estimate of 5i/m2 [|^. The idea is to apply an operator product expansion to 
the i?-meson matrix element of the time-ordered product of two flavour-changing 
heavy quark currents, and to equate the resulting theoretical expression to a 
phenomenological expression obtained by saturating the matrix element with 
physical intermediate states. This leads to sum rules, which can be used to 
derive inequalities for the B D^*^ transition form factors at zero recoil. In 
Ref. [0, such bounds have been obtained for the parameters ip and 5i/m2- They 



are 



?P > > , 
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Amcmb \3 



Ai + 2A^ 



2ml 



where 



pmm 
-P 



-(-A1-3A2), 



pmm 
-V 



The first relation in (^) implies that ||23| 

- Ai > 3A2 ~ 0.36 GeV^ 



(12) 



(13) 



excluding some of the values for the parameter Ai used in previous analyses 
of 5i/m2. It implies that the average heavy quark momentum inside the heavy 
meson is quite large, of order 600 MeV. Ball and Braun have calculated Ai using 
QCD sum rules and find -Ai = (0.5 ± 0.1) GeV^ [||]. Below we shall use 
Ai = —0.4 GeV^. We will comment on the (weak) dependence of our results on 
the value of Ai later. The upper bound for 6i/m2 in (p!l|) implies thatQ 77^^(1) < 
0.956. Of course, a crucial question is to what extent this inequality is saturated. 

make an "educated guess" that riA ^{l) = 0.89 ± 0.03 
= —(9.6 ± 3.0)%. However, the arguments presented to 



The authors of Ref. 
corresponding to Si/m'^ 
support this guess are not very rigorous. 

It seems more appealing to us to use the sum rules to constrain the hadronic 
parameters in (^. We first note that it is possible to derive two additional 
relations by interchanging pseudoscalar with vector meson states, corresponding 
to transitions of the type B* D^*\ Repeating the derivations of Ref. [0 for 
this case, we find 



> 2A2 - 0.24 GeV^ 



(14) 



5i 



I'm/ 



< 



( 1 



1 



/ /7mm 



nmm 



\2mc 2mf,J \2mc 2mf,J Arricmi, \3 



+ 



Ai + 2A2 



where ^1/^2 is obtained from 5i/m2 in (^) by interchanging £p with iy The first 
relation in ([1^) puts a bound on the parameter iy- To obtain further constraints. 



we use the fact that the above relations are valid for an arbitrary value of the 
mass ratio rric/mb. Comparing the second relation in ( pT]) with (^ in the limit 
TTib = rric, we find that 

Ag2>0. (15) 



^In Ref. this number is quoted as 0.94. 
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Figure 1: Allowed regions in the ip-iy plane for Ag'2 = 0.01 GeV^ 
(solid), 0.05 GeV^ (dashed), and 0.15 GeV^ (dash-dotted). We use Ai = 
-0.4 GeV^, in which case if"" = 0.02 GeV^ and i^'"" = 0.26 GeV^. 



Moreover, as long as rric < mb, it follows that 

777 . 777 

iv - iv" > — Up- C'") , 

rrif, — rric 

iv - if"" <—{ip- £g'") + Ag2 ■ (16) 

VTLc mi, — rric 

We are free to choose any value of the mass ratio mc/mf, between and 1 to make 
these relations as restrictive as possible. It is then straightforward to show that 

maxyip-if-- s/x^2; O} < < ^fi^^^W' + ^/>^^ ■ (17) 

For small values of Xq^, this relation implies a correlation between ip and iy, 
which is such that iy-ip^ i^'^'-if'' = 2A2 ^ 0.24 GeV^. This is illustrated in 
Fig. |l|, where we show the allowed regions in the ip-iy plane for different values 
of Ag'2. In total, we have thus identified three effects, which decrease Si/^'^ with 
respect to the estimate given in Ref. 0: a large value of (— Ai), a positive value 
of Ag2, and the fact that iy is likely to be larger than ip provided that Ag2 is 
small. 

To proceed, it is convenient to introduce new parameters 
i=l{(^v + ^p), 

s = lyv- ^r) + i^p - ^r)} = ^+li^i + A2) , (18) 
D = ^-{{iy- tr) - {ip - ^r)} = liiv- ip) - A2 , 

in terms of which 



The fact that S > imphes i > |(— Ai — A2) > A2. The inequahty (|T^) is 
equivalent to — -Dmax < D < -Dmax; where 

W^G^5-A|./4; 5>Ag./2. ^^^^ 

The main uncertainty in evaluating (|T9D comes from the unknown values of the pa- 
rameters i and Ag'2. As a guideline, one may employ the constituent quark model 
of Isgur et al. in which one uses non-relativistic harmonic oscillator wave 
functions for the ground-state heavy mesons, for instance ipsi'^) ~ exp(— i/icoT^), 
where /i = (1/mg + l/m;,)"^ is the reduced mass. One then obtains i = jm^ ~ 0.2 
GeV^, where we take ~ 0.5 GeV for the light constituent quark mass, cor- 
responding to the difference between the spin-averaged meson masses and the 
heavy quark masses. However, this estimate of i is probably somewhat too low. 
Lattice studies of heavy-light wave functions suggest an exponential behaviour 
of the form ipB{f) ~ exp(— K/xr) which leads to i = ~ 0.4 GeV^. We 
believe that values much larger than this are unlikely, since we use a rather large 
constituent quark mass nig. In fact, adopting the point of view that the sum 
rules for ip and iy are saturated to approximately 50% by the ground-state con- 
tribution |T0|, one would expect I ^ (— Ai — A2) — 0.28 GeV^, which seems a 



very reasonable value to us. In Fig. H, we show the allowed regions for 5i/m2 as 
a function of Ag'2 for two values of i. When — Ai is varied between 0.36 and 0.5 
GeV^, the resulting values for 6i/m'2 change by less than 1%. For all reasonable 
choices of parameters, the results are in the range —8% < Si/m'^ < —3%. Hence, 
we quote our new value as 

5i/^2 = -(5.5±2.5)%, (21) 

which is consistent with the previous estimates in Refs. |1^, ^ at the la 
level. A more precise determination of the parameter £ would help to reduce the 
uncertainty in this number. 

We conclude this section with a word of caution. Recently, it has been shown 



[P^ that the sum rules derived by Shifman et al. in Ref. |]10[ suffer from a renor- 
malon ambiguity; in other words, they do not obey the renormalization-group 
equation if the theory is regulated with a hard momentum cutoff. This is a se- 
rious problem, which has to be solved before these sum rules can be used with 
confidence in phenomenological applications. Here, we assume that the renor- 
malon problem can be cured without changing the form of the sum rules. 
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Figure 2: Allowed regions for 5i/m2 as a function of A^z for the two cases 
1= 0.2 GeV^ (solid) and 0.4 GeV^ (dashed). 



4 Prediction for the slope parameter 

In the extrapolation of the differential decay rate (|^) to zero recoil, the slope of 
the function ^iw) close to w = 1 plays an important role. One defines a slope 
parameter by 



l{nj)=a^){l-d\w-l) + 0[{w-lf]]. 



(22) 



It is important to distinguish from the corresponding slope parameter of 
the Isgur-Wise function. They differ by corrections that break the heavy quark 
symmetry. Whereas the slope of the Isgur-Wise function is a universal, mass- 
independent parameter, the slope of the physical form factor depends on loga- 
rithms and inverse powers of the heavy quark masses. On the other hand, 
is an observable quantity, while the value of depends on the renormalization 
scheme. 

To establish the relation between the two parameters, it is convenient to 
introduce in an intermediate step the axial vector form factor h^iiw) defined as 



{D*{vn*,e)\cj>'j,b \B{vb)) = {w + 1) HaAw) e*^ + . . . 



(23) 



where the ellipses represent terms proportional to f ^ or v^, . The relation between 
the physical form factors iiw) and hA^iw) is given in Ref. Defining a slope 
parameter in analogy to (^2|), we find 



Q 



6^ ^ 



1 



rriB 



3 mB — rriD* 



R 



2) , 



(24) 



where Ri and R2 denote certain ratios of the B D* decay form factors at zero 
recoil [0. In the heavy quark limit, Ri = R2 = 1, and the two slope parameters 
coincide. The symmetry-breaking corrections to these ratios have been analysed 
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in detail. Including both short-distance and I/ttlq corrections, which in this case 
can be calculated without much model dependence, one obtains Ri = 1.3 ± 0.1 
and R2 = 0.8 ± 0.1 ^. From (|2|), it then follows that Ap^ = -(0.22 ± 0.06). 
Recently, the form factor ratios -Ri and R2 have been measured by the CLEO 
collaboration, with the result that Ri = 1.30±0.39 and R2 = 0.64±0.29 @. This 
leads to Ag"^ = —(0.31 ±0.20), in nice agreement with our theoretical prediction. 

The next step is to relate the form factor hA^iw) to the Isgur-Wise function. 
The matrix element that defines the Isgur-Wise function in the heavy quark ef- 
fective theory is ultraviolet-divergent (for w ^ 1) and needs to be regularized by 
introducing a subtraction scale /i. To leading order in l/mq, the regularized func- 
tion ^(w,yu) is related to the physical form factor /i^^(w) by a renormalization- 
group-invariant Wilson coefficient function Cf, which contains the dependence on 
the heavy quark masses, and a universal function i^^^h containing the dependence 
on the renormalization scale JTU[ : 

hA^{w) = C^{mh,mc,w) Khi,{w, n) ^{w, fx) + 0{l/mQ) . (25) 

These functions are known to next-to-leading order in renormalization-group- 
improved perturbation theory. Using the explicit expression for i^hh('W^)A^) given 
in Ref. 0, we find that the physical slope parameter g\_^ is related to the slope 
parameter ^^(/x) of the regularized Isgur-Wise function by 

2 2/ N 16 8 /94 2\ "^(/^) 



d 

— In C^{mb,mc,w) 
aw 

^lnCf(m,,m„^ 
aw 



TT 

+ 0(l/mQ) 

UI = 1 

+ 0(l/mQ), (26) 

w=l 

,2 



where the last equation defines the /^-independent slope q of the renormalized 
Isgur-Wise function at next-to-leading order. Using the explicit expression for 
the Wilson coefficient given in Ref. [r^|, one finds that = + (0.21 ± 
0.02) -|- 0{l/mQ). An estimate of the l/mg corrections to this relation is model- 
dependent. We shall not attempt it, but instead increase the theoretical uncer- 
tainty to ±0.2. Hence, we obtain 

f = - (0.22 ± 0.06) - Q^± 0.2 . (27) 

Theoretical predictions for the renormalized slope parameter have been 
obtained from QCD sum rules, including a next-to-leading-order renormalization- 
group improvement. These calculations are tedious, since it is necessary to include 
two-loop radiative corrections to resolve the issue of scheme dependence. The 
complete calculation of these corrections has been performed in Ref. [^. It 
leads to = 0.7 ± 0.1 A similar result has been found by Bagan et al. [^. 
Based on (^), we thus predict 

^^ = 0.7 ±0.2. (28) 
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5 Summary 



The exclusive semileptonic decay mode B — >■ D*i v provides for the cleanest de- 
termination of the CKM matrix element Vcb- Heavy quark symmetry can be used 
to calculate the differential decay rate close to zero recoil in a model-independent 
way, up to small symmetry-breaking corrections, which can be analysed in a sys- 
tematic expansion in powers of as{mQ) and l/mg using the heavy quark effective 
theory. In this note, we have reconsidered and updated the analysis of these cor- 
rections. We find 7]a = 0.985 ±0.015 for the Wilson coefficient of the axial vector 
current, and 5i/m2 = —(5.5 ± 2.5)% for the power corrections to the normaliza- 
tion of the function ^{w) at zero recoil. The latter value is new and has been 
obtained by combining the existing approaches to estimate these corrections in a 
constructive way. Using these results, we predict 

=^^^(1) = 0.93 ±0.03 (29) 

for the normalization of the hadronic form factor T{w) at zero recoil. 

Three experiments have recently presented new measurements of the product 
I Vcb\ ^(1)- When rescaled using the new lifetime values r^o = (1.61 ± 0.08) ps 



and tb+ = (1.65 ± 0.07) ps ||3T|, the results obtained from a linear fit to the data 
are0 



0.0347 ±0.0019 ±0.0020 
Vcb\r]Ai{l) = { 0.0382 ±0.0044 ±0.0035 
0.0388 ±0.0043 ±0.0025 



CLEO [6], 
ALEPH 0, (30) 
ARGUS m. 



where the first error is statistical and the second systematic. Following the sug- 
gestion of Ref. P2|, we add 0.001 ± 0.001 to these values to account for the 



curvature of the function ^{w). Using then the theoretical result (p9D, we obtain 

I Vcb\ = 0.0399 ± 0.0026 (exp) ± 0.0013 (th) = 0.0399 ± 0.0029 , (31) 

which corresponds to a model- independent measurement of | Vcb\ with 7% accu- 
racy. This is by far the most accurate determination to date. 

We disagree with the conclusion of Ref. [0] that inclusive b ^ ciu decays 
would allow for a more reliable determination of | Vcb\- In this case, one has to 
make an assumption about the heavy quark masses that appear in the theoretical 
expression for the inclusive decay rate even at leading order. Moreover, it has 
been demonstrated that the extraction of | Vcb\ from inclusive decays suffers from a 
perturbative uncertainty of about 10%, due to unknown higher-order corrections 



in the expansion in as(mQ) Nevertheless, the most recent values obtained 



from the analysis of 6 — >■ c£i/ decays, which are 



I _ J 0.039 ± 0.001 (exp) ± 0.005 (th) ; measurements at T(4s), , . 
~ 1 0.042 ± 0.002 (exp) ± 0.005 (th) ; measurements at ^ ' 



^The ARGUS result has also been corrected for the new D branching fractions BS 
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are in excellent agreement with (|3T|). The theoretical uncertainty in these num- 
bers is larger than in the extraction from exclusive decays, however, and it is 
harder to control. 

Finally, we have related the physical slope parameter to the slope of the 
Isgur-Wise function and obtain the prediction = 0.7 ± 0.2 based on existing 
QCD sum rule calculations. It is consistent with the average value observed by 
experiments, which is g"^ = 0.87 ± 0.12 
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